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Abstract—Based on periodic solutions of the governing heat conduction equations in a single zone build¬ 
ing, computer software ADMIT has been developed for thermal simulation of buildings. Standard com¬ 
puter software, namely TRNSYS and SUNCODE, have also been used to simulate the same building 
under similar conditions. Simulations have been performed for three different climatic zones in India for 
light and heavy constructions under conditions of glazed/unglazed areas and ventilation rates. The results 
are presented in terms of the hourly variation of the room temperature. For insulated heavy construction, 
the results of different models are significantly different. This difference is due to the use of different 
approaches to solve the heat conduction equations. SUNCODE depends on the RC network approach 
and underestimates the heat losses. TRNSYS uses the transfer function approach, which is sensitive to 
the initially assumed value of the room temperature. ADMIT represents a quasi-steady-state periodic 
variation and is not suitable for transient variations. For insulated light buildings, the heat transfer 
mechanisms used in the mathematical models are not the governing factors. The models also differ in 
treating the penetration of solar radiation through a glazed window and the subsequent heat-transfer 
mechanism. For a south window and air changes in an insulated building, the results obtained by 
SUNCODE and ADMIT are in good agreement, but the results obtained by TRNSYS are considerably 
different. The reason for this needs detailed analysis. Copyright © 1996 Elsevier Science Ltd. 


1. INTRODUCTION 

Building simulation, from the point of view of 
calculating energy consumption for heating 
and/or cooling of buildings, has been of much 
interest in recent years. Task XIII of the 
International Energy Agency (IEA), with the 
aim of incorporating innovative materials, uses 
computer models to investigate the impact of 
advanced solar techniques on energy perfor¬ 
mance and thermal comfort in building design, 
components or systems in an advanced design 
(Eigdems, 1993). It was, however, found that 
most of these techniques cannot be directly 
modelled by the available simulation programs. 
The IEA task group, therefore, concentrated on 
a comparison of simulation models based on 
the calculation results of four typical cases and 
on composing information sets with regard to 
the modelling of the above-mentioned advanced 
technologies (Poel, 1993). The authors have 
been working with TRNSYS (1990) and 
SUNCODE (1985) for the past two years. It is 
found that modification of the existing codes 
needs a thorough understanding of several 
mathematical techniques, which is highly time 
consuming. In order to develop software where 
it is easy to incorporate several techniques of 
heating and/or cooling, the authors have used 
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the periodic solution technique of the governing 
heat conduction equations, resulting in closed- 
form solutions. The authors had the flexibility 
of incorporating many passive concepts, especi¬ 
ally for cooling, which are not yet available in 
the majority of the software. The results of the 
technique presented here are compared with the 
corresponding results obtained from TRNSYS 
and SUNCODE for light and heavy construc¬ 
tions and for different alternatives of physical, 
as well as climatic, parameters. 


2. MATHEMATICAL MODEL 


2.1. Periodic solution 


Heat transferred into a building can be 
assumed to be governed by the one-dimensional 
heat conduction equation: 


d 2 T(x,t) 8T(x, t) 

k ~e? — fC ~sT■ 


( 1 ) 


This has to be solved for various walls, roof 
and ground which are essentially flat geometries. 
Corresponding to Fig. 1, the boundary condi¬ 
tions, which have to be satisfied by the solution 
of eqn (1), are 


_ k mx 1 t) | 

(2a) 

dx | x=0 
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S„ = R 0 D„ + B„ (11) 


and 

Q n ^-R i E n + C n + R 1 R 0 D n -R 0 B n (12) 

where R a , R t are the surface film resistances, 
namely R 0 = \/h 0 and = 1 /h,. 

Solving the matrix eqn (10), the nth harmonic 
of the heat flux entering through a multilayered 
wall is written in terms of the temperatures on 
its two sides as: 


All time-varying functions, like solar radiation 
and ambient temperature, can be expressed in 
terms of a Fourier series, namely 

m = X /„ exp(inmr) (3) 


where to corresponds to the frequency of varia¬ 
tion, i.e. to = 27i/T p ; T p is the period of variation. 
For the nth harmonic of the temperatures T on 
and T r „ on the opposite sides of a composite 
wall, Pipes (1958) has shown that the temper¬ 
ature and energy cycles can be related to each 
other by the following matrix equation: 


eh; 
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where B in , C in and D in (i = 1 to m) correspond 
to the mth layer in the wall/roof and are given 
by: 

B m „ = cosh [(1 + i)a mn L m ] (5) 


^ sinh[(l +i)a m „L m ] 

C “-Ui + i)«™ <6) 

D m „ = -k m { 1 + i)ct m „ sinh[(l + i)a m „L m ] 

( 7 ) 
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Equation (4) can be written in the following 
form 


eh: -;]k a 

■C E1E1 


Qr n ~ 


S„ 7r„ — T on 

Qn 


(13) 


For the surface facing the ground, heat transfer 
is considered into a semi-infinite medium, R 0 is 
replaced by 


1 

k G (l + i) a G „ 


and T on = 0.0. 

In reality, the wall/roof and the floor receive 
solar radiation either on the external surface or 
on the internal (through a window) or on both. 
Its effect is incorporated into the above 
equations and the various heat gains through 
the corresponding walls/roof are obtained as 
given below: 

(1) The heat flux transmitted through the 
wall and roof is expressed as: 

Q = ^wl X 



Qn 


xexp (mart) + A wl X (Mb™) exp(incot) 
(14) 

where H WL is wall/roof area, / iw and /„ are 
radiation incident on inside and outside surfaces 
respectively, n denotes the nth harmonic, S„ and 
Q„ depend upon the thermophysical properties 
of walls or roof materials, h { and h 0 are the 
inside and outside convective heat-transfer 
coefficients respectively, aj and a c are absorptivi- 
ties of inner and outer surfaces respectively and 
oo is the angular frequency. 
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(2) The energy flux through the floor is 
obtained by: 

*- s -( T “ +!i fcf) 

Q = A F X ---exp(mcur) 

+ A f X (ot,/if„) exp(inojt) (15) 

where A F is the area of the floor and / if is the 
radiation incident on the floor. 

(3) Direct gain through the window is calcu¬ 
lated every hour and its effect on the indoor air 
temperature is calculated by assuming that 60% 
of the radiation is absorbed by the floor and 
8% is absorbed by each of the four walls and 
the ceiling. 

(4) The heat gain due to infiltration of air 
from the ambient into the room is calculated 
from the relation 

Q = C inf X (T an — T R „) exp(mcuf) (16) 

where C inf is the infiltration coefficient. 

(5) The expression for ventilation heat gain 
is given by: 

Q— X Cy(T an —T R „) exponent) 

= X E C vm (T a „- T RB )exp[i(n + m)cor] 

(17) 

where the ventilation term C v is assumed to be 
time dependent and is expressed in terms of the 
Fourier coefficient by: 

+ 00 

c v = X C vm exp (incur). 

(6) Heat transmission through the window 
glass and the door is expressed in terms of their 
respectives U values by: 

Q = A wn U X (-fan - T Rn ) exp (incot) 

(18) 

where A WN is the area of the window and U is 
the overall heat-transfer coefficient. 

(7) For the case when two rooms are adjacent 
to each other, heat transfer through the separat¬ 
ing wall (internal wall) is determined by the 
characteristics of the wall and is given by: 

Q = ^W L "f S " TRn Qn TR2n exp(incot). 

( 19 ) 


Room air temperature is essentially determined 
by the net amount of heat gain/loss by the room 
air through the building components. The same 
is given below, in the form of an energy balance 
equation given by: 

M t C p ^ X Tr n exp(/ncof)J = X Qj 

( 20 ) 

where M t is the thermal mass of room air, C p 
is the specific heat (J/kg K) and j corresponds 
to eqns (14)—(19). The choice about the use of 
number of harmonics usually depends upon the 
variation of the input parameters, i.e. ambient 
temperature and solar radiation. For smoothly 
varying inputs, even two harmonics, n = — 2 to 
+ 2, are sufficient. For fluctuating inputs six 
harmonics, n = — 6 to + 6, are found to be 
adequate. For n = 6 and n = 7, the results show 
a difference of ±0.001°C. Harmonics from n = 
— 6 to +6 are, therefore, considered and com¬ 
pared with coefficient of equal n. Equation (20) 
gives 13 equations which can be written in the 
form of a matrix equation given by: 

( 21 ) 

or 

iT Rn ^ 3xl =m; 3 \ l3 m l3xl . 

Solution of eqn (21) determines the different 
harmonic components of the room air temper¬ 
ature and these are combined to give the hourly 
variation of the room air temperature. 

TRNSYS and SUNCODE are essentially 
modelled with the same type of heat balances. 
The difference lies in the method of solving the 
heat conduction equations. In TRNSYS, the 
walls/roof are modelled according to the transfer 
function relationship of Mitlas and Arseneault 
(1971) evaluated numerically at equal time 
intervals. In SUNCODE the thermal response 
of walls is calculated by approximating the wall 
construction with a thermal netowrk. This net¬ 
work is solved using the explicit finite difference 
or Euler’s method. The program uses a central 
air temperature node for each zone defined by 
the user. The resulting air temperature is calcu¬ 
lated at each time step by rewriting the energy 
balance equations. 

3. RESULTS AND DISCUSSION 

The calculations have been performed for a 
single zone building of two sizes: 4 x 4 x 3 m, 
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and 6 x 4 x 3 m. The typical construction of 
roof and walls are given in Figs 2 and 3 for 
heavy and light construction, respectively; the 
thermophysical parameters used for the calcula- 

All dimensions are in mm 
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Fig. 2. Wall and roof details considered for a heavy 
construction single zone building. 



Fig. 3. Wall and roof details considered for a light 
construction single zone building. 


tions are given in Table 1. The calculations have 
been performed for a representative day of all 
12 months in three climatic zones of India 
(Bansal and Minke, 1988). The average climatic 
conditions of these three zones are given in 
Fig. 4. Additionally the following parametric 
variation is considered for the calculation: 

(1) no window, no insulation, no ventilation 
(NWNIO) 

(2) window, no insulation, no ventilation 
(WNIO) 

(3) window, no insulation, five air changes 
per hour (WNI5) 

(4) no window, insulation, five air changes 
per hour (NWI5) 

(5) window, insulation, five air changes 
per hour (WI5). 

For the case when a window is considered, it is 
assumed to be a part of the south wall and is 
10% of the wall area. When insulation of 5 cm 
polystyrene is considered it is assumed to be 
placed on the inner surface of the walls/roof, i.e. 
facing the room. For the sake of brevity the 
results have been presented for the months of 
January and June (typical winter and summer 
months, respectively) only. Figures 5-7 show 
the variation of room air temperature for 
different parametric variations for Delhi, 
Jodhpur and Leh, respectively, for a room size 
of 4 x 4 x 3 m. Figure 8 shows the same for 
Delhi for a room size of 6 x 4 x 3 m. 

Various effects can be seen in Fig. 5(a) and 
(b). For no window, no insulation and no air 
changes, hourly variation of the room air tem¬ 
perature calculated by (difference between the 
maximum and minimum temperatures) ADMIT 
shows the least fluctuations because this model 
is based on the quasi-steady-state periodic solu¬ 
tion approach. SUNCODE, which is based on 
RC network theory, estimates lower heat losses 
and therefore yields highest values for the room 
temperature. TRNSYS, based on the transfer 
function approach, is very sensitive to the ini¬ 
tially assumed value of the room temperature 
and also shows larger fluctuations than realized 
in practice. When a window is considered on 
the south wall, all the models yield higher 
temperatures, as expected, but the difference in 
the values of room temperatures between the 
three models is less than 2°C. Further consider¬ 
ation of five air changes per hour diminishes 
the difference still further. The difference 
between the calculated results using ADMIT, 
SUNCODE and TRNSYS never exceeds 1°C 
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Table 1. Thermophysical properties of different materials used for simulation 


S no. 

Material 

Density 

(kg/m- 3 ) 

Specific heat 
(J/kg K) 

Conductivity 
(W/m K) 

1 

brick 

1820.0 

880 

0.81 

2 

plaster 

1762.0 

840 

0.72 

3 

RCC 

2280.0 

880 

1.58 


mudphuska 

1622.0 

880 

0.52 

5 

broken tiles 

1820.0 

880 

0.81 

6 

insulation (expanded polystyrene) 

34.0 

1340 

0.035 

7 

soil 

1958.0 

840 

1.21 

8 

plywood 

530.0 

1300 

0.138 


except during the late hours, where TRNSYS 
still shows higher temperatures. While con¬ 
sidering a window and air changes in an 
insulated house, ADMIT and SUNCODE yield 
nearly the same room temperatures, whereas 
the temperatures obtained by TRNSYS are 
significantly lower. TRNSYS is therefore very 
sensitive to ventilation/infiltration of the ambi¬ 
ent air into the room. 

For lighter construction, as expected, the 
fluctuations in the room air temperature are 
higher. The difference in the room temperature 
values calculated by TRNSYS and SUNCODE, 
in this case, is marginal because the low capacity 
effect does not significantly affect the way heat 
transfer is calculated. For no window, insulation 
and five air changes, for example, the difference 
between the results is always lower than 0-5°C. 
As soon as the window is considered, the results 
of TRNSYS differ considerably (by about 3°C) 
with the corresponding values of SUNCODE 
and ADMIT. Obviously there is a difference in 
the manner in which the incident radiation is 
handled in various models. ADMIT calculates 
the incident solar radiation on windows, its 
penetration into the room, further reflection and 
absorption by other surfaces, etc., which actually 
happens. TRNSYS does take into account the 
absorption of reflected radiation but in a global 
fashion by weighing over the respective areas. 
This suggests that the penetration of solar radia¬ 
tion through a window and its distribution onto 
various surfaces need further investigations to 
modify the algorithm in TRNSYS and also 
perhaps in SUNCODE. 

Figure 6(a) and (b) shows the building’s per¬ 
formance for the hot and dry zone during the 
summer for heavy and light construction, 
respectively. For no insulation and no ventila¬ 
tion, the fluctuations in the room temperature 
are higher for lighter buildings as expected and 
the indoor air temperatures are well above the 
ambient temperatures. Consideration of five 


changes of air brings down the temperature. 
Providing an internal insulation and keeping 
five air changes helps to keep the room temper¬ 
ature around the average ambient temperature, 
as predicted by ADMIT and SUNCODE. 
TRNSYS shows higher fluctuations. Con¬ 
sidering a south window elevates the room 
temperatures, as predicted by ADMIT and 
SUNCODE, but the predictions by TRNSYS 
are still more optimistic, suggesting that 
TRNSYS needs to be investigated from the 
point of view of heat transfer and penetration 
of radiation through a window. Judkoff and 
Neymark (1995) corrected a bug in the transfer 
function (BID) module of TRNSYS causing 
insensitivity to thermal capacitance effects. In 
any case, all the models suggest that for hot 
and dry conditions, heavy construction with 
insulation on the external wall and avoiding a 
south glazed window is preferred for building 
construction. 

Figure 7(a) and (b) shows the hourly variation 
of room air temperatures for cold and sunny 
climatic conditions. The trends of variation pre¬ 
dicted by the models essentially remain the 
same. The overall physical behaviour is that the 
auxiliary heating and insulation of the buildings 
in these climatic conditions are necessary for 
thermal comfort. Lightweight construction with 
insulation behaves as good as a heavyweight 
building which is well insulated and possesses a 
window facing south. 

As the single zone size is increased, the 
variations in the results obtained by ADMIT, 
SUNCODE and TRNSYS are essentially the 
same (Fig. 8). 

4. CONCLUSIONS 

Different models use different mathematical 
techniques and therefore the results of building 
simulation from each of these models differ for 
different construction techniques and they are 
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Climatic Hot and dry Cold and sunny Composite 


Locations Jodhpur Leh New Delhi 

26°18'N 34°09'N 28°35'N 



Month 


Fig. 4. Climatic data for Jodhpur, Leh and New Delhi representation of hot and dry, cold and sunny, 
and composite climatic conditions. 


also sensitive to the climatic conditions. For indoor air temperature. All the three methods 
the case of a light insulated building with no differ significantly in dealing with the radiation 
window, the results obtained from different through a glazed window. Detailed comparative 
models are in good agreement, whereas for studies between simulation results and field 
heavy, multilayered walls/roof, the results are measurements are essential for further research 
significantly different. SUNCODE underesti- work. Moreover, incorporation of passive tech- 
mates the heat transfer (losses), while TRNSYS niques, especially for cooling purposes, should 
results are sensitive to the initially assumed be taken up for international use. 
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(A) 
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(B) 
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Fig. 7. Hourly variation of room temperature for Leh in January for (a) heavy construction (room size 
4 x 4 x 3 m), (b) light construction (room size 4 x 4 x 3 m). 


NOMENCLATURE h convective heat transfer coefficient (W/m 2 K) 

/ radiation flux (W/m 2 ) 

ach number of air changes per hour k thermal conductivity of the building material 

A area of fabric (m 2 ) (W/m K) 

c specific heat of the building material (J/kg K) L thickness of the elements (walls/roof) (m) 

Cm infiltration coefficient M, thermal mass of the room air (J/K) 

C v ventilation coefficient 








Temperature (°C) 
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P„, Q„, R„, S„ given by eqn (9) 

q rate of heat flow across any surface (W/m 2 ) 
Q heat flow across any fabric surface (W) 

R surface film resistance (m 2 K./W) 

T temperature ( C) 

71 time period (s) 
t time (s) 

V volume of room (m 3 ) 

X coefficient matrix of order 13 x 13 

Y column matrix of order 13x1 


Greek letters 

a absorptivity 

to angular frequency (per hour) 
p density of the building material (kg/m 3 ) 


Subscripts 

d door 
f floor 

i inside surface 
n nth harmonic 
o outside 
R room 

R 2 adjacent room 
v ventilation 
WL wall 
WN window 
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